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ABSTRACT 

Approximate analytical formulae are derived for the pulse profile produced by small hot spots 
on a rapidly rotating neutron star. Its Fourier amplitudes and phases are calculated. The pro- 
posed formalism takes into account gravitational bending of light, Doppler effect, anisotropy 
of emission, and time delays. Its accuracy is checked with exact numerical calculations. 
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1 INTRODUCTION 



During the past decade, coherent (or nearly coherent) oscillations 
in the ~ 200-600 Hz frequency range were discovered in the light 
curves of a number of neutron stars in low-mass X-ray bina- 
ries observed by the Rossi X-ray Timing Explorer (RXTE). In 13 
sour ces, these oscillations were discovered during X-ray bursts 
(see lStrohmaver & Bildster]|2006l . for a review), giving the name 
nuclear-powered millisecond pulsars to this class of objects. Seven 
transient sources, accretion-powered millisecond pulsars, showed 
coherent pulsations in the persist ent flux during t he outbursts last - 
ing a few weeks (see reviews by Poutanen 2006; Wiinands 1 20061) . 
In all these millisecond pulsars (hereafter MSP), the observed emis- 
sion is dominated by bright spots on the neutron star surface. Such 
bright spots are created either by a thermonuclear explosion ob- 
served as an X-ray burst or by an accretion flow channelled towards 
a magnetic pole. The pulse profile produced by the spot carries in- 
formation about the position of the spot, its size, spectrum and an- 
gular distribution of its emission and the gravitational field of the 
star. 

The X-ray pulse profiles observed from MSP are al- 
most sinusoidal. In SAX 11808 .4-3658 dGierlinski et~ai] 
120021; IPoutanen & Gieriinskil 120031) and XTE 11814-338 
dStrohmaver et alj |2003|) . there is a noticeable skewing of the 
profile which increases with photon energy. The skewness as 
well as the observed soft time lags probab ly result from Doppler 
boost of anisotropically em itted radiation (Gierlinski et al.ll2002l : 
Poutanen & Gierlinskil2003l) . Pulse profiles of SAX 11808.4-3658 
have been measured with high accuracy at diff erent photon ener- 
gies a nd well fitted by a theoretical model IPoutanen & Gierlinskil 
2003). This gave constraints on the neutron star radius 8 < R < 12 
km (assuming mass of 1.4-1. 6Mq) and on the inclination of the 
spin axis to the line of sight i > 60° . 

The poor photon statistics available for MSP normally does 
not allow one to study in detail the shapes of their light curves. 



Therefore, data analysis is often limited to the amplitude and phase 
of the fundamental Fourier harmonic, and sometimes higher har- 
monics may be analyzed. The Fourier technique helps use the 
available data to constrain the neutron star parameters (see e.g. 
i Miller & Lambll998l:lweinber g . Miller & Lambll200ll:lMuno et"al] 
120021 ; iGierlinski & Poutanen! 120051) . A number of recent papers 
were devoted to numerical calculations of the amplitude of pul- 
sations created by a spot (or two antipodal spots) for different radi- 
ation patterns, rotati onal velocities, inclinations, spot positions and 
their sizes (see e.g. Weinberg, Mille r & Larribl 1200 ll : | Mun o et al.1 
2002). Such calculations involve light-bending and Doppler effects, 
which made the problem complicated and required numerical cal- 
culations. Given the large number of parameters, the numerical ap- 
proach makes it difficult to understand the dependence of results on 
parameters and interpret the data. 

The purpose of the present paper is to develop an approximate 
a nalytical description of the problem. Using the simple formalism 
of lBeloborodovl j2002L hereafter B02) for light bending, we derive 
analytical formulae for oscillation amplitudes which demonstrate 
how the observables depend on parameters of the problem. 

The plan of the paper is as follows. In Section [2] we intro- 
duce our notations and summarize the exact numerical method of 
the light-curve calculation. In Section[3]we discuss the approximate 
description of light bending, introduce an approximate formula for 
time-delay effects, describe possible classes of MSP, and finally de- 
rive approximate analytical formulae for the pulse profiles and the 
corresponding amplitudes and phases of the Fourier series. Using 
our analytical formalism we compare the approximate Fourier am- 
plitudes to the exact results and investigate the effects of anisotropy 
of the emission pattern and the Doppler effect on the pulse profile 
in Section|4] 



2 PULSE FROM A SPOT ON A SPINNING STAR 



* E-mail: juri.poutanen@oulu.fi (JP); amb@phys.columbia.edu (AMB) 



Consider a small spot on the star surface. Its area measured in the 
corotating frame is dS', and its instantaneous position in the fixed 
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Figure 1. Geometry of the problem. Dotted curve shows the photon trajec- 
tory. 



lab frame is described by the unit vector n that points to the spot 
from the star center (see Fig.[TJ. The angle between n and the line 
of sight is denoted by ip. We are interested in photons emitted by 
the spot that propagate along our line of sight at large distances 
from the star (where gravitational bending becomes negligible). We 
denote the unit vector along the line of sight by fc, so that 



cost/; — k ■ n. 

As the star rotates, k ■ n varies periodically, 

cos ip — cos i cos 9 + sin i sin 6 cos <f>, 



(1) 



(2) 



where i is the inclination angle of the spin axis to the line of sight, 
6 is the spot colatitude and <j> = 2nvt is the rotational phase of the 
pulsar; v = P" 1 is the pulsar frequency, and t = is chosen when 
the spot is closest to the observer. 

Angle ip measures the apparent inclination of the spot to the 
line of sight, which is different from the true inclination because 
of the light bending effect. We denote the initial direction of the 
emitted photon by feo and the true emission angle by a, so that 



cos a = feo ■ n. 



(3) 



Emission angle in the corotating frame is denoted by a'. It 
differs from a because of relativistic aberration (see derivation in 
the Appendix) 



cos a' = 5 cos a, 

where 5 = 1/7(1 — P cos£) is the Doppler factor. Here 7 
/3 2 ) -1//2 and (3 = v/c is the spot velocity, 



2nR 



sin 6 — p ccl sin f 



(4) 



(5) 



/3 eq is the velocity at the equator and £ is the angle between the 
spot velocity and feo. Here u = rs/R, rs = 2GM/c? is the 
Schwarzschild radius; M and 7? are mass and radius of the star. 
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Figure 2. (a) The bolometric black body flux as a function of the observed 
phase. Dotted curve is for a slowly rotating star ignoring all relativistic ef- 
fects. The pulse profile where gravitational light bending is accounted for 
is shown by the dashed curve. Dot-dashed curve gives the profile modified 
by the Doppler boost and aberration for a neutron star rotational frequency 
v = 600 Hz. Solid curve accounts also for the time delay. Thin solid curve 
is a pulse profile produced using the approximate Fourier amplitudes de- 
rived in Sect. |3 (b) Doppler factor 8 as a function of the observed phase for 
v = 600 Hz. We take i = 9 = 45°, R = 2.5r s , M = 1.4M Q in this 
example. 



The pulsar frequency has been corrected for the redshift yl — it. 
One can show that £ is related to a, %[>, i and 4> by (see Appendix), 



cos£ = 



sin Tp 



(6) 



For power-law spectra (observed for example in accretion- 
powered MSP), we assume that the energy and angular dependen- 
cies of the spectrum emitted by the spot may be separated as 



I E '{a) = io(l + hcosa')E'- {T ~ 1 \ 



(7) 



where h does not depend on photon energy E. If the power-law 
spectrum is produced by thermal Comptonization by electrons of 
temperature T c , this condition would be satisfied if the maximum 
Doppler shift AS ~ 2-kvR/c sini sin (9 is smaller than the 
typical relative energy change in a single scattering AE/E ~ 
AkT e /m e c 2 , which translates to (v/60 Hz) sini sin 8 < 
fcT c /16 keV dViironen & Poutanenl [2004b . Even for X-ray burst 
spectra one expects that h varies slowly with energy, so that equa- 
tion ifTJ still can be used. 

The observed spectral flux at a distance D from the star is then 
given by (see derivation in Appendix), 



F E = (l-u) T/2 8 T+3 l' E (a')> 



dcosa dS" 



d cos ip D 2 



(8) 



Expression for the bolometric flux may be obtained as a special 
case of equation l[8j by setting T = 2, 



F = (1 - u) <f I' (a) 



d cos a dS" 



(9) 



d cos ip D 2 

These equations take into account the special relativistic effects 
(Doppler boost, relativistic aberration) as well as general relativistic 
effects (gravitational redshift and light bending in Schwarzschild 
geometry). 

For further analysis we use pulse profiles normalized to Fo = 

Jo£-< r - 1 )(l_ tt )C r + ! 0/ a dS'/D 2 : 
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Figure 3. Light bending in Schwarzschild metric. The solid curves are the 
exact results using equation )Alt and the dotted lines is the approximation 



F((j)) = 8 l+i {l + hS cos a) cos a 



d COS Q 



(10) 



1 — u d cos i/j 

The flux is zero if cos a < 0. For the antipodal spot, we substitute 
f? — > 7r — f3 and <f> — > n + 4>. 

Expression d 1 0b does not account for time delays resulting 
from different paths travelled by photons emitted at different phases 
4>. The delays become significant only for very fast-rotating pulsars. 
In Schwarzschild metric the maximum time delay for a neutron star 
of M — 1.4M Q is At ~ 7 x 1CT 2 ms (almost independent of 
compactness of the star M/ R). This gives at most a 5 per cent cor- 
rection to the arrival phase for a rotational period P — 1.5 ms. The 
flux at observed phase </> b s is F((j>obs) = F((/> b B — A<f>) with 
phase delay Acj> = 2nvAt computed using (AT7\ and jAT8l from 
the Appendix. The effect of the photon arrival time contraction (or 
stretching) on the observed flux is already accounted for by one of 
the Doppler factors. The effects of gravitational bending, Doppler 
boost, and time delays on the pulse profile are shown in Fig. [2] 

Fourier series of the pulse profile is given by 



where 

A ~- 

A n = 
B n - 



A + / ][A„ cos(n^ obs ) + B n sin(n0obs 



(11) 



1 

2^ 
1 

W Jo 
1 



F((/> ob s) d(f> 



F((f>obs) cos(n</> obs ) d^obs, n > 1 



F((f>obs) sm(n(fi ohB ) d</>obs, n > 1. 



(12) 



An alternative form of Fourier series is written in terms of ampli- 
tudes c n and phases A(f> n 



F(^obs) = / ) Gn cos[n(^ obs + A<j> n )], 



where 



C n = s/Al + Bl, taa(nA^ n ) = -B n /A n . 
The exact Fourier series are calculated numerically. 



(13) 



(14) 



Figure 4. Time delays in Schwarzschild metric. The solid curves are the 
exact results using equation jA17t for R = 2, 2.5 and 3rg (from top to 
bottom). The dashed curves show the results of an approximate formula 
I18K and the dotted line is the approximation H9\ . 



3 ANALYTICAL APPROXIMATION 
3.1 Light bending and time delay 

The expressions for the observed flux {§} and ((9) can be signifi- 
cantly simplified if one uses analytical formula for light bending 
and time delays. This will allow us to also obtain analytical expres- 
sions for the pulsation amplitude and Fourier harmonics in Section 

E31 

iBeloborodovl d2002h showed that the relation 

cos a ~ u + (1 — u) cosi/i (15) 

descr i bes light bending with high accuracy (see also Zavl in et"al] 
1 19951 ; iLeahv &THll995l , for discussion of approximations). The 
accuracy of equation l |15t is shown in Fig. [3] (see also B02). For a 
star with R = 2rs, the accuracy is better than 10 per cent, while for 
R = 3rs the error does not exceed 3 per cent. The spot is visible 
to the observer when 

cos-)/' > cosVWx = — w/(l — it). (16) 

The exact expression for the time delays l |A17t may also be 
approximated by a simple formula. Expanding it in Taylor series 
with 1 — cos a as a small parameter and using also the expansion 

2 



COS Oi I U 2 



1-u \ 112" 

(obtained from eq. [2] in B02), where y 

uy 



(17) 



cos ip, we get 



At 



y{l + ^[l + y(l/3-u/14)]} R/c 



Keeping only the leading term, we find 

At = y R/c. 



(18) 



(19) 



This approximation is better than 20 per cent accurate for most 
emission angles and compactnesses (Fig. 0- Only in the extreme 
case of R — 2rs for large bending angles (and large delays), the 
error increases to 35 per cent. Since the time delays themselves pro- 
duce a small effect, equation l |19l > is sufficient in most calculations. 
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Figure 5. Classes of visibility for hot spots on neutron star surface at the inclination i- spot colatitude 6 plane. The classes are determined by the relative 
positions of /x m ; n = cos(i + 6), /imax = cos(i — $), cos i/> ma x and re = | cos ^>max|- Left panel: classes for one spot. In class A defined by /U m i n > — re the 
spot is always visible. In class B defined by £t m i n < — re < /i max the spot is not visible for a fraction of period. In class C defined by ^ ma x < — K the spot is 
never visible. Right panel: classification for two antipodal spots as given in B02. In class I, fj, mlIL > re, the primary spot is visible all the time and the antipodal 
spot is never seen. In class II, —re < fj, mln < re < /i mM , the primary spot is seen all the time and the antipodal spot also appears for some time. In class III, 
Mmin < — re, the primary spot is not visible for a fraction of period (and then only the antipodal spot is seen). In class IV, —re < fj, mln , (i maI < re, both spots 
are seen at any time. Dotted and solid lines correspond to the exact (f mM = 129? 4 and approximate cos t/i max = —uf (1 — u) (giving t/i max = 131?8), for 
u = 0.4 (i.e. R = 2.5r s ). 



3.2 Pulsar visibility classes 

Analytical light bending formula U5\ allows one to introduce a 
simple classification of the light curves according to the relative 
positions of /ti m i n = cos(i + 9), /i max = cos(i — 9), cos^ max 
and ft = | cos^ ma x|. For a single spot, three classes exist (see left 
panel Fig.[5J. In class A, defined by fimin > — ft, the spot is always 
visible. When ft < — /x max (class C), the spot is always invisible. 
For /i m in < — k < /i max (class B), the spot is visible during the 
where 



pulsar phases \(p\ < 

COS0p = (cOS^mas 

and we defined 



cos i cos 9) I sin i sin ( 



U = (1 — u) sini sin^, 
Q = u + (1 — u) cos i cos ( 



-Q/U, (20) 



(21) 



Pulsars with two antipodal spots are divided into four classes 
(shown in the right panel of Fig. [5] see also B02). In class I, corre- 
sponding to Umin > k only the primary spot is visible all the time. 
In class II, — ft < /i m in < ft < /x max , the primary spot is always 
visible, while the antipodal secondary spot appears during phases 
(f> s < (f> < 2tt — (j> s , where 



cos 0s = 



COS 'i/'max + COS 1 COS 9 2u — Q 

sin i sin 9 U 



(22) 



The primary spot disappears for a fraction of the period in class III, 
/imin < — ft, and then only the antipodal spot is seen. And finally, 
in class IV (—ft < /Lt m i n , Atmax < ft), both spots are seen all the 
time. 



3.3 Fourier series 

Our aim is to obtain simple analytical expressions for the Fourier 
amplitudes and phases characterizing the pulse profile. We can sim- 
plify expression dlOt by using approximation to the light bending 
formula J 15b : 



F(<p) 



rr+3 



( 1 + h8 cos a) cos a, 



where now 

cos a = Q + U cos ( 



(24) 



Doppler factor 8 depends on /9 cq . In the leading order of /3 cq <t; 1 
this dependence is given by 



o ~ 1 — p c „ sm i sin 8 sin < 

smifi 



(25) 



We further approximate sina/sini/> ~ y/1 — u which becomes 
exact ata< 1 (cf. the cosine relation 1151 ). This gives 



5 « 1 - T sin ( 
where 



T = /9 cq \/l — u sini sinf? <JC 1. 



(26) 



(27) 



Hereafter, we keep linear terms in T and neglect higher order terms. 
Then, substituting 8 into expression l |23t , we obtain the flux (as a 
function of (f>): 



F{4>) = a o + ^2[ a n cos(n0) + b n sin(n0)], 

n=l 

with non-zero coefficients 

-i2 . rr 2 



ao 

Ol 

bi 

a z 
b 2 

b 3 



Q + h(Q^ + Er/2), 
{l + 2hQ)U, 

- Q(3 + r) + /i'- >J 



+ 



u 2 



(4 + r) 



T, 



hU 2 /2, 

-[(l + 2ftQ)(4 + r) 
- A -±LhTU 2 . 



TV 



(28) 

(29) 
(30) 

(3D 

(32) 

(33) 

(34) 



(23) 



Similarly to equation l !14b . we can define Fourier coefficients c n 
and corresponding phase lags. Coefficients l|29|l-l|34t are good ap- 
proximations to the exact Fourier coefficients A n , B n of the flux 
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F(cj> b s ) as long as the spot is visible at all phases (through the en- 
tire period of rotation). Alternatively, the pulse profile can be writ- 
ten as a cosine series 

F(<t>) = Q + h{Q 2 + U 2 /2) 
U(l + 2hQ) 



COS fl 

hU 2 /2 



cos(<j) + Ci) 



(35) 



- f- cos[2(<£ + £,)] + ^-hTU 2 cos[3(<?(> + tt/6)], 

cos 2C,2 4 



where 



. , T (3 + r)Q + (4 + r)h(Q 2 + U 2 /4) 

tanCl ~ u TT2KQ ' (36) 



tan 2(2 



T (4 + T)(l + 2ftQ)-l 



V h 
The ratio of the harmonic to the fundamental 



(-2 



hU 



cos (1 



ci 2(1 + 2hQ) cos2C 2 



(37) 



(38) 



grows with the anisotropy parameter h and sin i sin 9. 

If the spot disappears from the visibility zone for part of the 
period (classes B, II, and III), the Fourier series becomes infinite. 
The Fourier coefficients of the light curve can be computed from 
the coefficients j29t-<34t: 



a = — y^a k s k (<j)p), 

TV ^— ' 

fc = 

a'n = — / ak[s„-k((f>p) + Sn + fe(0p) 
TV — ' 

b' n — — y bk[s n -k{4>p) — S„+k(<j>p) 

TV • ■* 



(39) 



where so(c/>) = 1 and s n (<j>) — sm(n<f>) / (ruf)) . 

Including the time delays further modify the expansions. We 
are interested in the Fourier coefficients of the function F((f> b s ) = 
F(4»obs — A0(0)). We can calculate the phase delays relative to 
the photons arriving from the star element closest to the observer 
(with impact parameter b = 0) using equation H9\ : 



Acj>((j>) fa A0(0 obs ) w i> - Tcos0 obs . 



(40) 



where $ = /3 eq \/l — w(l — cos i cos 6 1 ) = (1 — Q)/3 cq /\/l — M - 
Keeping only the first term in Taylor expansion, we arrive at 



cosn(0 — A0) = cos n(j> + nA(f> sin n<f>, 
sinn(0 — A0) = sin ruj) — nAcj) cos n<f>. 

The Fourier amplitudes for F(<f) b B ) are found as follows 



a n = a 



(41) 



(42) 



T 



b„ + $na„ — ^ [(n — l)a n -i + (n + l)o n +i] • 



where we neglected products T6 n ,$fe n oc T 2 . If the spot is in- 
visible for a part of the period, one should use coefficients a' n , b' n 
instead of a n ,b n . For the antipodal spot, we substitute 9 — > n — 8 
and (j> — > 7T + 4>. Thus we still can use the expressions d29t-ll34ll. 
where Q is replaced by Q B = u — (1 — u) cos i cos 6 and the sign 
of the odd terms 01, 61, 63 is changed. 



4 RESULTS 

We have checked the accuracy of the analytical formalism 
(Sect. 13.31 ) by direct comparison with the exact numerical calcu- 
lation (Sect.O. The accuracy depends mainly on the compactness 
of the star u = rs / R and its rotational frequency v. Part of the er- 
ror comes from the light-bending approximation dl5t ; this error in- 
creases for stars with large u (Fig. [3}. Besides, we made the approx- 
imation sin aj sin tp fa (1 — u) 1 ^ 2 in equation J25b for the Doppler 
factor 8. This introduces an additional error which becomes notice- 
able for fast rotators. Errors generally grow at higher v because we 
neglected the terms quadratic in /3 eq in all formulae. 

As an example, we show in Figures [2] and [6] the results for 
a fast and compact rotator with v — 600 Hz and u — 0.4 
(R/rs = 2.5). In this case, the pulse profile reconstructed from 
the analytical Fourier amplitudes is very close to the exact profile 
(Fig.[2](. The accuracy of analytical approximation for amplitudes 
Co and Ci is better than 3 per cent, while C2 is 15 per cent ac- 
curate (see Fig- [6j. The phases fi and £2 (eqs. 1361 and 1371 . and 
the phases obtained from eqs. 1391 and 1421 ) are accurate within 
0.2 rad. For slower rotation, e.g. v — 300 Hz, the error of ana- 
lytical approximation decreases to ~ 5 per cent. Only amplitude 
C3 (which is smaller than Ci and C2 and more sensitive to the 
neglected terms in the analytical expansion) has a significant error, 
our formulae underestimate C3 by about 50 (20) per cent for ^=600 
(300) Hz. 

Even for an extremely compact star with u = 0.5, the analyt- 
ical C\ and C2 have a good accuracy: they are typically 10-20 per 
cent smaller than the exact values. Only in the cases with one spot 
at large colatitudes 6 > 120° and small inclinations i < 60°, the 
amplitudes are underestimated by a factor of 2, because of extreme 
gravitational bending. 

Using the analytical formalism one can understand the be- 
haviour of Fourier amplitud es and their ratio ob t ained previously b y 
numerical calculations ^Weinberg et"alll200ll ; iMuno et alj|2002l) . 
Below we investigate separately the effects of anisotropy and fast 
rotation. 



4.1 Effects of anisotropy 

Let us assume a slowly rotating star (i.e. take 5 = 1 and o b s = <f>) 
and investigate the effects of anisotropy h 7^ of the source 
emission 1(a) — Io(l + hcosa) on the pulse profile. For 
example, radiation from an optically thin source (slab) is more 
beamed along the stellar surf a ce and described by h < (see e.g. 
IPoutanen & Gierlinskil 120031 : Iviironen & Poutanenl |2004|) . while 
h ~ 2 characterizes radiation from an optically thick electron- 
scattering dominated atmosphere. 

As will be seen from equations below, anisotropy (i) intro- 
duces harmonic 122 oc h, (ii) changes the global structure of the 
profile (e.g. two maxima may appear and cj> = may become a 
minimum), and (iii) leaves the profile symmetric. 

We will use the following expression for the flux from a single 
(or primary) spot, 



Fp{4>) = cos a p [1 + h cos a p ] , 
where 



(43) 



cos Qp = u + (1 — it) cos ip = Q + U cos 4>. (44) 
If secondary (antipodal) spot is present, the corresponding flux is 

F E (4>) — cos a B [1 + h cos a B ] , (45) 
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Figure 6. Fourier amplitudes of the pulse profile as functions of the spot colatitude for three inclinations i = 30° , 60° , 90° . Left panel: one spot. Right panel: 
two antipodal spots. Solid curves show the exact numerical results and dashed curves - the analytical approximation. Inclinations and corresponding visibility 
classes are marked above the curves. Parameters are M = 1.4Mq, R = 2.5rg, r = 2, v = 600 Hz, h = —0.7. 



where (using cos ip s = — cos tjf) 

cos as = u — (1 — u) cos ip — 2u — Q — U cos 4>. (46) 
4.1.1 One spot 

First, let us consider the simple case of black body spots (i.e. h = 
0). The pulse profile from a single spot then has almost exactly sine 
shape (B02) 



F p ((f>) = cos a p = Q + U cos 4>. 



(47) 



This immediately gives the non-zero Fourier cofficients ao = Q 
and ai = U for pulsars in class A. In class B, the Fourier coef- 
ficients can be computed using equations ([39} . The peak-to-peak; 
amplitude A = (F max - i ? min)/(i ? max + Fmin) is given by: 



{U/Q, class A, 
1, class B. 
0, class C. 



(48) 



The peak-to-peak amplitude coincides with the ratio |ai|/ao and is 
\pl larger than sometimes quoted rms amplitude if the light curve 
is a pure sine profile. 

The pulse profile is given by 

/ rr2 \ tj2 

F{4>) = Q+h ( Q 2 + — J +U(l+2hQ) cos <f>+h— cos 20.(49) 



One sees that anisotropy introduces no phase shift in the harmonics, 
and the pulse profile remains symmetric about <j> — and <f> = 
7T, however the Fourier amplitudes are changed. The amplitude of 
the fundamental ci increases with positive h and decreases if h 
is negative. The fundamental can even completely disappear when 
h = — 1/(2Q). It is proportional to U and therefore ci behaves 
close to sin i sin 6 (see also Fig. [6] left panel for class A and right 
panel for class I). 

Anisotropy / introduces harmonic cos 2cf> in the pulse. 
The amplitude C2 of this harmonic can be quite large if h is large 
(in absolute value) and is proportional to sin 2 i sin 2 6 (see eqs. 1321 , 
1331 and lower panels in Fig.[6). The ratio of amplitudes C2 and c\ 
is given by 



£2 _ 02 _ hU/2 

ci ~ en ~ 1 + 2hQ 



oc ft sin i sin 9. 



(50) 



If h ^ —1/2, the pulse profile has a maximum at (j> — and 
a minimum at (f> = n. The corresponding pulse amplitude is 



.4 



U(l + 2hQ) 
' Q + h(Q 2 + U 2 )' 

If h < — 1/2 and the conditions 

U +l/2/i 



cos(i + o) < r\ = — - 



1 - u 



< cos(i 



(51) 



(52) 
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are satisfied, then the pulse has minima at both <j> = and (/> = ir. 
Two maxima 



4/i 

appear at phases corresponding to 
1 + 2hQ 



COS <p = COS </>!=—- 



2hU 



The global minimum of the pulse profile is at <f> — n. 



-fmin ( 



tt) = (Q-U)\l + h(Q-U)]. 



(53) 



(54) 



(55) 



If h < —1/2 and cos(i + 8) > rj, the pulse minimum is at 
<f> — and its maximum at = n. When cos(i — 9) < r\, then the 
minimum is at <j) = and the maximum is at <j> = 7r. 

In class B, the phase 4> — t is not visible. The condition for 
additional maxima d52b is the same as in class A and these maxima 
are always in the visibility zone, cos <f>i > cos cj> p . 

The above analysis is also applicable to pulsars with two an- 
tipodal spots in class I, when only one spot is visible. Because 
class I is smaller on the 9-i plane than class A, the existence of 
the maxima at <f> — \cf>i\ (eq. 1541 ) requires a stronger condition on 
h: -l/4u <h< -1/2. 



4.1.2 Two antipodal spots 

First consider black-body spots (h = 0). In class IV, the pulse pro- 
file is flat with ao = 2u, i.e. there are no pulsations. In class II, the 
pulse profile consists of a single sinusoidal pulse and a flat plateau 
appearing when the antipodal spot also becomes visible (see eq. 
1221 ). In class III, there are two pulses with two plateaus in be- 
tween (symmetric about phases cj> — and n). The primary is seen 
at \(f>\ < (j>p, while the secondary appears at cos(/> < cos(/> s . The 
corresponding pulse amplitudes A for all these classes are given by 
(B02) 



A = I 



U/Q, 
U + Q 



2u 



U + Q + 2u 
0, 



cos(i — 6) — u 
cos(i — 6) + u ' 



class I, 

classes II, III, 
class IV. 



(56) 



The anisotropy h 7^ modifies significantly the pulse profiles. 
The analysis is simplest for class IV, where two spots are always 
visible. Then the total flux is 

F{<j>) = F p (4>) + F s (4>) =2u + h[2(Q-u) 2 + 2u 2 + U 2 

+ 4(Q -u)U cos <j) + U 2 cos 20] (57) 
and the Fourier amplitudes are 

2h(Q-u)U 



£1 _ Oi 

c ao 

£2 _ 02 

Co ao 



u + h[v? + (Q-u) 2 + U 2 /2Y 

hU 2 /2 

u + h[u 2 + (Q -u) 2 + U 2 /2]' 



(58) 
(59) 



The dependencies ci/co oc h sin 2i sin 29 and C2/C0 oc 
h sin 2 i sin 2 9 are found not only in class IV, but also in classes II 
and III (right panel, Fig.[6]l. Note that the importance of the second 
harmonic 



£2 _ U 
ci 4(Q - u) 



■ tan i tan 9 



(60) 



is independent of h and grows with i and 9. 

If i + 9 sj 90°, the pulse profile has two extrema at (j> — and 
(f> — tv. The pulse peak-to-peak amplitude is 



A = 



2(Q-u)U\h\ 



, + h\u 2 + (0 - u) 2 + U 2 } ' 



(61) 



If i + 9 > 90°, there are two additional extrema at cos<j> — 
cos tf>2 = — (Q — u)/U = — cot i cot 9 with 



F(fa) =2u(l + hu). 



(62) 



These extrema are minima if ft > and maxima if h < 0. Here 
<f> = is a global maximum for h > and a global minimum for 
h < 0. The amplitude of the pulse with four extrema grows with 
\h\. 



A = 



(1-it) 2 cos 2 (i-(9) \h\ 



(63) 



2u + h[2u 2 + (1 - u) 2 cos 2 (i - 9)] ' 
It is maximum at the boundary of class IV when cos(i — 6) = k. 

4.2 Effects of fast rotation 

The effects of fast rotation are illustrated in this section with a sim- 
ple case of black body spots (h — 0, no anisotropy). Furthermore, 
we assume that the spots are always visible, i.e. consider classes A, 
I, IV. 



4.2.1 One spot, classes A and I 
The flux from one spot is given by 

Fp(» = 5 3+r cosa p (64) 
= Q + Ucoscj>- (3 + r)OTsin0- UTsm2<j> 

= Q+ cos(0 + Ci) + ^-±-^l7Tcos[2(0 + jr/4)], 

COS C]_ ^ 

which only slightly deviates from a sinusoidal shape. The phase 
shift 



tanO = (3 + T)Q 



T 



u yr 



:(3 + T)Q 



(65) 



is of the order of /3 oq 1. The main dependences of the Fourier 
amplitudes on i and 8 are similar to the anisotropic case considered 
in Section 14.1.11 with ci behaving close to sin i sin 9 and C2 oc 
sin 2 i sin 2 9 (compare eqs. 1301 and 1311 as well as eqs. 1321 and 
1331 ). The ratio of harmonics is 



£2 _ (3 + ryr 

ci ~ 2 



(3 + r) 



it sin 1 sin ( 



(66) 



We see that the Doppler effect produces the first harmonic and in- 
troduces a phase shift between the fundamental and the harmonic, 
which skews the profile. The phase shift is proportional to /3 eq - 

We note, however, that the ratio c-2/ci depends linearly on 
/3 eq (which is a small number) and sin i sin 9. We conclude that 
for a single spot C2/C1 <§C 1, i.e. the Doppler effect alone cannot 
introduce a strong additional harmonic to the signal (unless V is 
large). By contrast, an anisotropic source can make it easily (see 
eq. 1501 ). In combination with anisotropy, the Doppler effect sig- 
nificantly modifies the pulse profile. 

The time delays slightly reduce the phase lag to tan £1 = 
/3c q [(4 + T)Q — — u an d produce a third harmonic with 

a small amplitude 03 oc UT 2 oc /3 2 q . 
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4.2.2 Two spots, class IV 

The contribution of the secondary (antipodal) spot to the flux is 
given by 



Fs(4>) 



c3+r 



cos a s = 2u 



(3 + r)(2it-Q)Tsin< 



U cos 4> 

3 + r 



UT sin 2(j), (67) 



where S s ~ 1 + T sin (f>. The total flux is given by 

F(» = F p (4>) + F.(4>) 

= 2u~ (3 + F)T[2(Q -u)sm(j> + Usm2(j)]. (68) 

The time delays introduce the effects of the order /3 2 q which are 
ignored here. 

The Doppler effect modulates the flux, so that the profile is 
not a plateau anymore. The Fourier amplitudes ci and c-2 grow lin- 
early with the rotational frequency and the spectral index V. The 
amplitude of the fundamental, 



ci 3 + T 



(Q — u) T (X /3 cq sin 2i sin 29 



(69) 



c u 

is maximum at the boundary of class IV where cos(z — ff) = k, 
while the amplitude of the harmonic, 

c 2 3 + r , 



- UT oc /3 eq sin 2 i sin 2 l 



(70) 



Co 2u 

is maximum when i = 6 — 90°. The dependences on i and 9 are 
similar to the anisotropic case from Section 14.1.21 This explains 
the fact that when both anisotropy and fast rotation are present, the 
behaviour remains the same (see right panels in Fig. [6}. The ratio 
of the harmonics is 

~ = ^77T" — T = \t&nitan9. (71) 
ci 2(Q -u) 2 

It differs by a factor of two from the anisotropic case <60l l. and 
becomes large at large inclination i and spot co-latitude 9. 

If i + 9 ^ 90°, the Doppler-boosted pulse from two antipodal 
spots has two extrema at 



cos — cos ( 



V8 + cot 2 i cot 2 9 — cot i cot 9 



>0, 



(72) 



with maximum at ~4>+ an d minimum at <J) + . If i + 6 > 90°, there 
are two additional extrema at 



v8 + cot 2 i cot 2 9 + cot i cot 9 „ ,_„. 
cos<p = cos0_ = < 0, (73) 

with <f>- being the maximum and —cj>- being the minimum. The 
global maximum and minimum remain at =F0+. 



5 SUMMARY 

We have derived an analytical approximation for the pulse profile 
produced by small spots on a neutron star surface, its Fourier ampli- 
tudes and phases. The exact profile is reproduced with good accu- 
racy by our formulae even in the case of very fast rotation (e.g. with 
the error of ~ 15 per cent for rotational frequency v = 600 Hz). 
For slower rotation of 300 Hz, the accuracy improves to ~ 5 per 
cent. The main advantage of the analytical formalism is that it 
shows the dependence of the pulse profile and its Fourier series 
on the parameters of the pulsar. 

The proposed formalism can be used to obtain constraints on 
the neutron star parameters and position of the hot spot from the 
amplitudes of oscillations observed during X-ray bursts in some 



low-mass X-ray binaries, as well as from the X-ray pulse profiles 
of the accretion-powered and rotation-powered pulsars. Our results 
can be further extended to study the energy dependence of the pro- 
file as well as the effects of the finite spot size. 
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APPENDIX A: EXACT CALCULATION OF OBSERVED 
FLUX 

Al Light bending and Lorentz transformations 

The exact relation between a and rj) in Schwarzschi ld geometry 
(i.e. light bending) is given by (e.g. lMisner et alll973l) 



b 2 r' 2 ( r ) 



1 -1/2 



where b is impact parameter, 



(Al) 



(A2) 



u = rs/R, rs = 2GM/c? is Schwarzschild radius; M and R 
are mass and radius of the star. The maximum bending angle i/Wx 
corresponds to a = n/2. The visibility of the spot is defined by 
a condition cos a > 0, or alternatively by ip < t(i m „ = %p(a = 

7T/2). 
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As a photon emitted at angle a with respect to the spot nor- 
mal n propagates to infinity, its direction changes from fco near the 
stellar surface to fc at infinity, so that cos a = feo • n changes to 
cos ip — k ■ n. The relation between fco and k may be written as 



feo = [sin a k + sin(^) — a) n]/ sin ip. 



(A3) 



At any moment of time, we can introduce an instantaneous non- 
rotating frame x,y,z with the j/-axis along the direction of the spot 
motion, x-axis along the meridian towards the equator, and z-axis 
along the normal to the spot. In this static frame, 



feo 



(sin i cos 9 cos <j> — cos i sin 6) , cos £, cos a ) , (A4) 



where 

(3 sin a /3 sin a ... cx 
cos £ = — ■ fco = ■ k — sm i sin <p . (A5) 

p sini/i p shif/; 

In the frame comoving with the spot (with y-axis along the spot 
motion, z-axis along the local normal), the unit vector along the 
photon momentum is obtained from the Lorentz transformation: 

((sin i cos 8 cos <f> — cos i sin 9) sin a / sin ip \ 
7(cos£-/3) , (A6) 

cos a / 



where 7 = 1/ ^/l — f3 2 and the Doppler factor 

7(1 — /3 cos £) 
Using equation 1 IA6I 1, we obtain 

cos a = 6 cos a. 



(A7) 



(A8) 



where E/E' — 5\/l — u. Here I' E , (a') is the intensity computed 
in the frame comoving with the spot. For the bolometric intensity, 
one gets 



(SVT^) 4 I' {a). 



(A13) 



If the radiation spectrum can be represented by a power-law 
I' E i{a') oc E'~^ r ~ 1 ^ with a photon spectral index Y which does 
not depend on the angle a' then 



l' E ,{a') = l' E { a '){8Vl~^) 1 



(A14) 



This approximation is equivalent to the assumption of a 
weak energy dependence of the angular distribution (see 
IViironen & Poutanedl2004h . 

The observed spectral flux (eq. |A9t now reads 



dF_E = (1 — u) 1 ' 2 5 i I E i(a) cosa 



d cos a dS' 



(A 15) 



dcos-;/) D 2 ' 

where we have used the aberration formula dA8b . Substituting equa- 
tion l lA14t . we recover equation {8j. 
The bolometric flux is given by: 



dF = (1 — u) 5 i' {ex) cos a 



dcosa dS" 
dcosi/' D 2 



(A 16) 



Thus, the flux from a rapidly rot ating star differs by a factor 5 5 from 
that from a slowly rotating star jPoutanen & Gierl ihski 2003). Two 
powers of S come from the solid angle transformation, one from the 
energy, one from the photon arrival time contraction, and the fifth 
from the change in the projected area due to aberration. Aberration 
may also change the specific intensity since it has to be computed 
for angle a' in the comoving frame. 



A2 Observed flux 

The observed flux from the spot at photon energy E is 
dF E = Ie dfi, 



(A9) 



where Ie is the specific intensity of radiation at infinity and dil is 
the solid angle occupied by spot with area dS' on the observer's 
sky. The solid angle can be expressed in terms of the impact pa- 
rameter 



dQ = bdb dtp/D 2 



(A10) 



where D is the distance to the source and ip is the azimuthal an- 
gle corresponding to rotation around line of sight (vector k). The 
impact parameter b depends on ijj only, but not on tp. 

Using equation dA2t and the facts that dS = R 2 d cos ip dip 
and dS" cos a' — dS cos a (since the spot area projected on to 
the plane perpendicular to the photon propagation direction, i.e. a 
photon beam cross-section, is Lorentz invariant), one gets 



dQ = 



dS" cos a' 1 d cos a 



D 2 



Uud cos ip 



(All) 



In the limit of weak gravity » < 1, this gives the usual formula 
dO = dS'cosa'/D 2 . 

The combined effect of the gravitational redshift and Doppler 
effect results in the following relation betw een the monochro- 
matic observed and local intensities (see e.g. iMisneretalJl 1973k 
iRvbicki & LightmarJI 19791) : 



(A 12) 



A3 Time delays 

Finally, we write down here the formula describing the time delays. 
The delay is caused by different travel times of emitted photons to 
the observer, depending on the position of the emitting spot. A pho- 
ton following the trajector y with an impact para meter b is lagging 
the photon with b = by ( Pechenick e t all 19831) : 



cAt{b) 



poo 

Jr t 



dr 



r s /r 



1 



('-?) 



-1/2 



1 } .(A17) 



For a given pulsar phase <f>, we compute angle ip, then we find the 
corresponding emitted a and the impact parameter using formu- 
lae JAlb and JA2b . and compute the corresponding delays At(b) 
with equation JA17b . We then construct a one-to-one correspon- 
dence between the pulsar phase <p and the photon arrival phase to 
the observer <p \, s = <f> + A(f>, with the phase delays 



A0(0) = 2nuAt[b(4>)] 



(A 18) 



For analytical work we can also use the relation <j> = o b s 

A(j}(<j}) ~ 0obs — A0(<^ o ba)- 



